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^H Abstract 

^^ We discuss the dependence of J/V' production on the charged particle multiplicity in proton-proton 

rvj collisions at LHC energies. We show that, in the framework of parton saturation or string interaction 

^^ models, the hard J/i/) production exhibits a significant growth with the multiplicity, which is stronger than 

f^ linear in the high density domain. This departure from linearity, that should affect any hard observable, 

rr^ applies for high multiplicity proton-proton collisions in the central rapidity region and is a consequence of 

tn the parton saturation or the strong interaction among colour ropes that take place at LHC energies. Our 

(-H assumption, the existence of coherence effects present in proton-proton collisions at high energy, can also be 

checked by studying the particular shape of the probability distribution associated to the J /xp production. 
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I RHIC [U l^l [21 13] and LHC [HI E] data on heavy-ion collisions have shown several important features 

^^ which indicate the formation of a high density partonic medium with characteristic properties as the low 

y^ shear viscosity and high opacity. Since the energy density achieved in high multiplicity events produced in pp 

;_( collisions at LHC is comparable to the reached density in CuCu central collisions or AuAu peripheral collisions 

at ^smn = 200 GeV, it is pertinent to wonder about the possibility to obtain a similar high density medium 

which would be reflected in experimental observables, similar to heavy-ion collisions. An illustration of this 

would be the predicted ridge structure [8j [9] observed by CMS collaboration [10] in pp collisions. Also, the 

eventuality that other observables, as long range rapidity correlations [HI [121 US]; energy loss [14], or the 

elliptic flow [151 116] . are measurable in pp collisions, has been reckoned with in different frameworks. 

We address here to the J /^j production in high multiplicity pp collisions. We will show that the rise of 
J /ip production in the highest multiplicity events observed by the ALICE collaboration [17] can be naturally 
explained as a consequence of string interaction or parton saturation. This feature is particularly important, 
since, due to the absence of nuclear effects in this case, it can be use as a baseline in order to disentangle different 
mechanism that are expected in heavy-ion collisions, as the J/ip suppression due to sequential dissociation [TS] 
or the J/tp enhancement due to the recombination of uncorrelated c and c quarks [H] . 



Our main assumption, shared by many of the mentioned approaches, is the fact that, in high-energy hadronic 
coUisions, all projectiles, be it protons or nuclei, have finite spatial extension and thus collide generically at 
finite impact parameter by means of elementary parton-parton collisions. We may consider the color ropes or 
flux tubes - strings - as the fundamental variables of our description. They are formed in each parton-parton 
collision and they constitute the elementary sources of particle production. In this string framework, the number 
of parton-parton collisions is reflected as the number of produced strings, Ng. These strings have non- negligible 
transverse size, of the order of 0.2 -;- 0.3 fm and different space-time rapidities, and they can interact - overlap -, 
so reducing the effective number of sources, in particular in which concerns soft particle production. 

Consider now the hard particle production. The number of initially produced J/ip, "j/i>i can be taken as 
proportional to the number of collisions, in analogy to any hard process. In the string-like models, this number 
corresponds to the number of produced strings, Ng. On the other hand, the rapidity multiplicity distribution 
dN /drj -mainly soft- is not proportional to the number of collisions, but mostly to the number of participants. 
This reduction can be considered as a consequence of shadowing (20] , parton saturation |21j or string interactions 
- percolation - [25 . 

In the string percolation approach, the multiplicity distribution is given by 

^ = F(p)7V,Mi (1) 

drj 

where jii corresponds to the multiplicity of a single string in the rapidity range of interest, Ng is the number of 
produced strings and F{p) corresponds to the damping factor induced by the string interaction. 



I - e^P 

np)-^^^. (2) 

Note that, within the damping factor, 1 — e^^ represents the fraction of the total area that is covered by 
strings. The interaction among strings and the consequent damping factor is a function of the the string 
density, p — °'^° , where cro is the transverse area of one string, ctq — ^rrp, rg ^ 0.25 fm, and a corresponds to 
the transverse area of the collision. So the existence of many strings -directly related to the number of available 
partons- effectively screens the charged particle multiplicities, which agrees qualitatively with the concept of 
saturation. 

On the other hand, assuming the proportionality between the number of produced J/ij) and the number of 
collisions, 

(3) 



ij/V _ ^s 



< nji^ > <Ng> 

it is possible to obtain the relation between the charged particle multiplicity and the number of produced J/')p, 
that, accordingly to eqs. ([2| and ([s]), will obey 



dN , .1/2 

A) _ / '^J/i/J 



<^> V<"j/^> 
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where < p >=< Ns> — . (4) 

a 



dr\ 

At low multiplicities, where the number of strings < A^^ > is small, the above equation gives rise to the 
linear dependence 



dN 



^ r, , ~-^ ^ dN ^ ■ ^^> 

< nj/v > <^> 



On the contrary, at high niuhiphcities, the bracket in the right hand side of eq. Q can be approximated by 
< p >^^/^. One obtains in this case 

C dN \ 2 
^^ ■ (6) 

Thus the Unear dependence obtained previously for \ow multiphcitics, eq. ([5|, changes to an squared dependence 
when high muhiphcity events are at play. 

In order to compare with the available pp experimental data, we take a = cr^^^^i = 70 mb [2 3) as the transverse 
area of the collision. The number of strings, < Ng >, can be obtained from the SFM code [Mj, a Monte Carlo 
code based on the quark gluon string model, similar to the Dual Parton Model [25] or EPOS [26]. Moreover, 
the value oi Ng can also be calculated analytically [27j, < Ns >= b+ {2 — b) (-p) , where A = 0.2-^0.3, b = 1.37 
and the low energy threshold ^/st = 10 GeV. We obtained, for the central rapidity region, < Ng >— 16, while 
in the forward rapidity region the number of strings is smaller, < Ng >= 8. The reason for this difference is 
the fact that, while the long strings -stretched between valence quarks and diquarks of the colliding protons- 
cover most of the rapidity range thus contributing to both central and forward rapidity production, the short 
strings -stretched between sea quarks and antiquarks- are mostly created in the central rapidity region only. 

In Fig. 1 we show our results for both the central and forward rapidity range, together with the experimental 
data from the ALICE Collaboration |l7j. We observe a good agreement. The main uncertainty of our evaluation 
concerns the value of < iV^ >, which can induce an uncertainty in the < p > value of 20%, of the same order 
of the last experimental point error. Moreover, it would be most interesting to measure ^^'"^ > ^^ larger 

dN 

multiplicity, ^^ ~ 8 -f- 10 to check the behaviour of data. Note that for this high pp multiplicities, the 
values of the parton densities p are comparable to those obtained in CuCu central collisions or AuAu peripheral 
collisions at RHIC energies, where J/V' melting is observed. If experimental data were significantly lower than 
our curve, this would be an indication of J/i/; suppression in absence of any nuclear effect. 

Our assumption, the existence of coherence effects present in pp collisions at high energy, can also be checked 
by studying the probability distribution associated to the J/'0 production, Pc{n). This distribution is given by 

Pe{n) = -^^^P{n) , (7) 

where P{n) corresponds to the minimum bias probability distribution. The above equation refers to the well- 
established universal relation between minimum bias distribution P and the multiplicity distribution P^, associ- 
ated to the production of a rare event c [IHl HH] ■ This formula is valid in any model of particle production based 
on the superposition of the contributions from elementary partonic interactions -as it is the case in most of the 
multiple scattering models-, when the kind of events trigerred is shadowed only by events of the same kind, and 
not for the whole of events. In addition to this requirement, these events must have a small cross section. In 
other words, one can say that the rare events are those produced only by one of the elementary interactions, the 
probability of it occurring twice in one collision being negligible. This is the case if one is triggering on a heavy 
particle like the J/ip, and it would also be the case when studying T production or multiplicity distributions 
associated to W, Z^ production in pp collisions. The validity of eq. (|7|) has been checked for several cases in 



A 




' 


>> 




,' 


5 10 


_ 


1 


> 




/ 


z 


- 


/ 


■a 




/ 


_v 


- 


1 


>. 




' 


■a 




/ 


■-> 




,' 


Z^ 8 


1 


1^ 


■a 


/ 




6 


^ 


K- /^ 


4 


- / /^ 




2 

n 


— 1 — 1 — 1 — — 1 — 1 — 1 — — 1 — 1 — 1 — — 1 — 1 — 1 — L 


1 1 1 1 1 1 L 



4 5 6 

dN,,/dTl/<dN,,/dTi> 



Figure 1: Our results for the central {\y\ < 0.9, dashed line) and forward (2.5 < y < 4, dotted line) rapidity 
range, together with the experimental data for the central (circles) and forward (squares) rapidity regions 
from the ALICE Collaboration j/7,. The linear behaviour (solid line) is also plotted. 



AA and pp collisions [28j[29]. Moreover, a different trend for the multiplicity distributions for interactions with 
and without charm, in agreement with our expectations, was already obtained more than twenty years ago by 
NA27 Collaboration [50] in pp collisions. This difference was interpreted to be due to the more central nature 
of collisions leading to charm production. This is also the case in our approach, when taking the shape given 
by eq. 0. 

From the experimental minimum bias distribution P{n) measured by ALICE [17 for the rapidity range 
\ri\ < 1.0, and using the relation established in eqs. p| and mK for the computation of '^^^^'^ as a function of 
the multiplicity, one can calculate the probability distribution for the J/V", Pc{n), accordingly to eq. (l7|. Our 
results for the central rapidity region are shown in Fig. 2, compared to the minimum bias distribution P{n) 
in this region. A particular shape, characterized by a reduction of the probability for low multiplicities and an 
increase at high multiplicities when compared to the minimum bias distribution is obtained. Moreover, we have 
checked that, when using PYTHIA 6.4 in the Perugia 2011 tune [3T1 El] for the computation of ^„''^'''y as a 
function of the multiplicity, accordingly to the result presented in [171 \22\ in the central rapidity region, and 



applying ([t]), the opposite behaviour resuh^ 
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Figure 2: Minimum bias experimental distribution P{n) (solid line) compared to our obtained Pc{n) distribution 
for the J/-0 (dashed line) and to the Pc{n) distribution obtained when using PYTHIA for the J/tp production 
(dotted line). 



In conclusion, we have reproduced here the rise of J/ip production in the highest multiplicity events observed 
by the ALICE collaboration in pp collisions. This increase, more pronounced in the central rapidity region, may 
be a consequence of the formation of a high density medium in pp collisions at LHC energies. In this case, the 
linear dependence of J/ip production on the charged particle multiplicity obtained for low multiplicities -where 
the parton densities are smaller-, changes to an squared dependence when high multiplicity events are at play, 
due to the high string densities. 

This behaviour can be checked by studying the probability distribution associated to the J/^p production, 
compared to the minimum bias probability distribution. A particular shape, characterized by a reduction of 
the probability for low multiplicities and an increase at high multiplicities when compared to the minimum bias 



* According to 1171 . the PYTHIA result exhibits a decrease of the J/ip multiplicity with respect to the event multiplicity, contrary 
to data. 



distribution would be obtained. 
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